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We have detected an x-ray nebula around the binary milhsecond pulsar B1957+20. 
A narrow tail, corresponding to the shocked pulsar wind, is seen interior to the 
known Ha bow shock and proves the long-held assumption that the rotational 
energy of millisecond pulsars is dissipated through relativistic winds. Unre- 
solved x-ray emission likely represents the shock where the winds of the pul- 
sar and its companion collide. This emission indicates that the efficiency with 
which relativistic particles are accelerated in the post shock flow is similar to 
that for young pulsars, despite the shock proximity and much weaker surface 
magnetic field of this millisecond pulsar. 



Millisecond pulsars are old neutron stars (typically ~ 3 Gyr) that have been spun up to a 
rapid rotation rate (< 25 ms) by accretion of material from a binary companion (Tll^. After the 
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accretion phase, they appear as radio pulsars with surface magnetic field strengths of ~ 10^ G, 
which, combined with their older ages and rapid rotation rates, means that they form a separate 
population from younger pulsars. 

PSR B 1957+20 is the second fastest-spinning pulsar known, with a rotation period (T) of 
1.6ms and rotational spin-down luminosity {4} E = 1 x 10^^ ergs s^^ The pulsar is in a 9.16- 
hour binary orbit with a low-mass companion star. The wind of the companion star eclipses 
the radio emission for ~10% of every orbit. The PSR B 1957-1-20 binary system provides an 
excellent opportunity to study the wind of a weakly magnetized, recycled neutron star. The 
wind is ablating, and may eventually evaporate, the low-mass companion star. Ablation and 
heating of the companion star Q are believed to be caused by x- or 7-rays generated in an 
intrabinary shock between the pulsar wind and that of the companion star. Meanwhile, the high 
space velocity of the pulsar (> 220 km s^^) fES, as it moves through the interstellar medium 
(ISM), generates sufficient ram-pressure to confine the pulsar wind and results in the formation 
of a bow shock. Upstream of this bow shock Ha emission is generated where neutral material is 
swept up and coUisionally excited; this Ha nebula absorbs 1-10% of the spin-down energy Q. 
A reverse or termination shock decelerates the pulsar wind and is located between the pulsar 
and the bow shock. 

The intrabinary and termination shocks both provide diagnostics of the pulsar wind. Ob- 
servations of the interaction of young pulsar winds and the ISM have shown that they power 
x-ray emitting synchrotron nebulae which are typified by the Crab nebula ©. High-resolution 
observations of such pulsar-wind nebulae have shown that the composition of the wind, char- 
acterized by the ratio of the Foynting flux to particle energy flux a, is kinetic energy dominated 
(cr ^ 0.003). However, for the evolutionarily distinct millisecond pulsar population, the com- 
position, efficiency, and geometry of the pulsar wind remain unknown. ROSAT observations of 
PSR B 1957-1-20 provided the only previous constraint, and they have been interpreted as indi- 
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eating that its wind is different from that of the Crab pulsar ©. However, these data lacked the 
spatial resolution and sensitivity to determine the cause of the emission. 

We have undertaken a 43kilosecond dT^ observation of the millisecond pulsar B 1957+20 
using the Chandra X-ray Observatory. The brightest x-ray source in the field (Fig. 1) is coinci- 
dent with the pulsar position. A tail of x-ray emission is seen extending from the pulsar to the 
northeast by at least 16" with a position angle opposite to the pulsar's proper motion direction 
of 212° ®. 

Striking confirmation of the association with the pulsar and the ISM shock comes by com- 
paring this x-ray tail with the Ha bow-shock nebula (Fig. 2). In bow-shock nebulae, material 
at the termination shock is swept back by the ram pressure and forms a cylindrical tube aligned 
with the proper motion direction and interior to the bow shock il3\\14li . The morphology, di- 
rection and location of the x-ray nebula (Fig. 2) therefore indicate that it corresponds to the 
termination shock and, in combination with the enclosing Ha emission, demonstrates the ex- 
pected double-shock nature of the pulsar's interaction with the ambient medium. 

Young pulsar winds are thought to be relativistic, and correspondingly generate non-thermal 
emission by synchrotron and/or inverse Compton processes (I75li<5l) . However, nothing is known 
about millisecond pulsar winds except that the scale of the observed Ha nebulae are at least 
consistent with their having relativistic winds (I77t . A different wind might be expected because 
of the reduction and possible alteration of the nature of the magnetic field during the accretion 
phase. We must therefore also consider the possibility that they drive slower winds like that of 
the Sun, resulting in shock-heated thermal x-rays. If the spin-down luminosity of a millisecond 
pulsar is carried away as kinetic energy in an outflow then E = Q.5MV^, where M is the mass 
loss rate and is the wind speed. Because the pulsar remains and is at least 1 Gyr old we know 
that M < 10^^ g s^^, and thus Vw > 10^ cm s^^. The temperature that we measure {kT ^ 1 keV, 
where k is the Boltzmann constant) by fitting a spectrum expected for shock-heated gas to the 
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energy distribution of the 82 counts recorded in the tail region is inconsistent with such high 
velocities HM . 

The emission must therefore be from some form of non thermal process, either synchrotron 
or inverse Compton emission, as seen in wind nebulae around young pulsars. In either case, 
the x-ray emission requires a population of relativistic particles in the pulsar wind. Thus, the 
detection of a distinct x-ray tail provides direct evidence that millisecond pulsars lose their 
rotational energy through relativistic winds. 

It is highly unlikely that diffuse shock acceleration dTPt is the mechanism that generates 
the emitting particle population in this nebula, because it is difficult to accelerate particles in a 
relativistic flow through this mechanism (l20b . As in other pulsar wind nebulae, the acceleration 
presumably takes place through some other mechanism possibly involving the role of heavy 
ions 

Because the pulsar and the companion are separated by only 1.5 x 10^^ cm, the intrabinary 
shock, formed where the pulsar and companion winds interact, will be located in a strong mag- 
netic field (much stronger than that at the shock in the Crab nebula). This intrabinary shock is 
therefore a potential source of unresolved synchrotron emission at the location of the pulsar (|22b 
and can be used to determine the wind characteristics at the shock front. The x-ray luminosity 
in the shock is dependent on both the post shock magnetic field strength and a (HI. We therefore 
consider two possibilities for the wind composition of PSR B 1957-1-20: either dominated by ki- 
netic energy (a = 0.003 as seen in the Crab nebula) or magnetically dominated (a ^ 1). We 
can describe the process by which the spin-down energy of the pulsar is converted into x-ray 
emission at the intrabinary shock by {9): 

ffjAeL^ = fradf-yfgeomE (1) 

where is the spectral x-ray luminosity in a band of width Ae ~ 1 keV centered on a pho- 
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ton energy e = 1 keV. The geometric factor fgeom defines how much of the wind interacts 
with the companion, is the fraction of the intercepted spin-down energy flux that goes into 
accelerating electrons with Lorentz factor 7 corresponding to e = 1 keV, and frad is the ra- 
diative efficiency of the corresponding synchrotron emission. The fraction of the unresolved 
x-ray emission that is due to synchrotron emission produced at the intrabinary shock is /;,, the 
remainder presumably being emission produced by the pulsar itself. 

Table 1: Properties of the PSR B 1957-1-20 system ("S"), and comparison with the Crab nebula 
and its central pulsar The fraction of the wind of PSR B 1957-1-20 that interacts with the 
companion star wind is fgeom ~ ''^e/(4a^) ~ 0.02, where a is the orbital separation. This fraction 
is a lower limit, because the pulsar wind is most likely focused into the equatorial plane (I24II25I) 
which is also probably the orbital plane of the binary system The companion star would 
therefore intercept more of the spin-down energy from the pulsar than if it has a spherical wind. 



Parameter 


B 1957-1-20 


Crab 




(a = 0.003) (a > 1) 


Nebula / pulsar 


Spin period (ms) 


1.6 


33.5 


Surface magnetic field (10^ G) 


1.4 


3.8 X 10^ 


Distance (kpc) 


1.5 


2 


Age (yr) 


> 2 X 10*^ 


948 


Distance to shock (cm) 


- 1.5 X lO^i 


~ 3 X 10^^ 


E (10^5 erg s-^) 


1.0 


4.4 X 10^ 


h 


~ 0.5 




fgeom 


> 0.02 




B(G) 


2 12 


^ 10"^ 


tflow (s) 


5 1.7 


~ 3 X 10^° 


trad (s) 


700 45 


~ 5 X 10" 


frad 


0.007 0.04 


0.06 


£ (7 10^) 


1 keV 


70 fim 


L, (7 ^ 10^) (erg s'l keV"i) 


2.6 X 10^0 


5.8 X 10^° 


A (7 ~ 10') 


< 0.09 < 0.02 


0.04 



The flow immediately downstream of the intrabinary shock is expected to undergo Doppler 
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boosting as it passes around the companion (E^ . Thus, one expects the x-ray emission at orbital 
phases before and after eclipse to be enhanced by up to a factor 2.2 depending on the flow speed 
and the degree of absorption and/or scattering by the companion wind (E^ . In contrast, the x- 
ray emission at the eclipse (orbital phase 0.25) may be reduced because of obscuration of the 
shock by the companion star. If we can measure these variations then we can determine /;,. The 
lowest and highest count rates in the folded light curve (Fig 3.) are during and immediately 
after eclipse, respectively, with comparatively low probability of this variation being by chance 
(see caption of Fig 3.). If this apparent orbital modulation is genuine and corresponds to a 
modulation of the x-ray flux from the intrabinary shock by a factor of 2.2, then we find that 
fb ~ 0.5. 

The post shock magnetic field strengths ® corresponding to our two limiting values of a 
(0.003 and ^ 1) are listed in Tabled the corresponding Lorentz factors of the synchrotron 
emitting relativistic electrons for both cases are 7 = 2.4 x 10^ (e / BY^"^ ^ 10^ at e = IkeV, 
where B is the post shock magnetic field strength. For each case the emitting region is taken to 
be the radius of the radio-eclipse region ^ Ve = 5 x 10^° cm; for a = 0.003 the flow speed is 
Vfiow = c/3, whereas for cr ^ 1 we expect vjiow = c (0). The corresponding residence times 
are given by tfiow = vjiow/re (Table [B, and the radiative lifetimes of the emitting electrons 
are Uad = 5.1 x 10^/(7!?^). We can thus compute the radiative efficiency (l26b in each case 
as frad = (1 + Uad/tfiowY^- With /g^om > 0.02 (Table [B and ~ 0.5, we can infer that 
< 0.09 (a = 0.003) or < 0.02 (a > 1) for Lorentz factor 7 ~ 10^ 
We directly compared the properties of the shocked wind of PSR B 1957-1-20 with those 
derived for the Crab nebula with the same Lorentz factor and find that 0.04 for electrons 
with 7 = 10^ in the Crab, which, regardless of the assumed value of a, is similar to that which 
we derived for PSR B 1957-1-20. The available evidence therefore suggests that despite being 
subject to a prolonged evolutionary process that has altered its magnetic field by many orders 
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of magnitude and having a shock that occurs much closer to the pulsar, this millisecond pulsar 
generates a wind for which the efficiency of relativistic particle acceleration in this post shock 
flow is similar to that seen for the winds of young pulsars. New Chandra data on this and other 
pulsars (I27II2^|EPI) are providing the first detailed observational input into studies of relativistic 
flows, particle acceleration, and magnetohydrodynamic shocks. 
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Figure 1: The full resolution Chandra image from the ACIS-S3 detector in the energy 
range 0.3-10.0 keV. The source at (J2000) RA 19^59°'36!75±0!01; Dec +20°48'15'.'0±a'l 
is coincident with the proper motion-corrected radio timing position of PSR B 1957+20, RA 
19^59™36!75788±O:00005; Dec +20°48'14'.'8482±a'0006 (6). Four of the next brightest x- 
ray sources in the field are circled and the source labeled 1 is coincident to better than O'.'l 
with the Tycho II star 1628-01136-1 {11). In an aperture of radius 1'.'5 centered on the pulsar 
position we detect a background corrected total of 370 ± 20 counts from the pulsar in the en- 
ergy range 0.3-10.0 keV. Binning the data in energy such that each bin contains a minimum of 
30 counts, we fit an absorbed power law resulting in a best fit with hydrogen column density 
Nh = (1.8 ± 0.7) X 10^^ cm"^, photon index F = 1.9 ± 0.5, and an unabsorbed flux density 
Fx,c = 6 X 10~^^ ergs s'^ cm^^ (0.5-7.0 keV) corresponding to an isotropic x-ray luminosity 
ofL^^c = ^T^D'^F^^c = 1-6 X lO^^Di.s ergs s^^ (0.5-7.0 keV), where D1.5 = D/(1.5kpc) is 
the pulsar distance as derived from its dispersion measure (17^ . Counts in the tail region were 
summed in a 16" x 6" box enclosing the tail and aligned with the proper motion direction. After 
background correction, we detect a total of 82 ± 9 counts from the tail region. Again fitting an 
absorbed power law and assuming that the Nh and F are similar to the values above we derive 
an unabsorbed flux density F^^^t = 9 x 10~^^ ergs s"^ cm^^ (0.5-7.0 keV). 
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Figure 2: The same Chandra data as shown in Figure 1, smoothed to a resolution of 5" (con- 
tours) and overlaid on an Ha image obtained from Taurus Tunable Filter service mode obser- 
vations on the Anglo Australian Telescope on 2000 Aug 3. The images were aligned to within 
~0^'2 through use of USNO 2.0 stars. The x-ray tail is located well inside the boundaries of the 
Ha emission and also close to its symmetry axis. The x-ray contour levels are shown at 0.9, 
1.2, 5.3, 35.0 and 78.8% of the peak x-ray surface brightness. The optical residuals correspond 
to incompletely subtracted stars. No optical counterpart was found to the x-ray source associ- 
ated with the contours located in the northeastern comer of the field on a Digitized Sky Survey 
image. 
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Figure 3: The light curve of the x-ray emission (0.3-10.0 keV) from the point source associated 
with PSR B 1957+20 folded through use of the known orbital ephemeris (D. J. Nice, private 
communication, 2001). The errorbars correspond to one-sigma Poissonian errors. Radio eclipse 
occurs at orbital phase 0.25. The mean count rate, excluding phases in the range 0.15-0.35 
where variations in the x-ray flux are expected, is (8.4±0.8) ctsksec^^. Orbital phase bins 0.25 
and 0.35 each correspond to a total observing time of 6600s, and thus we expect to measure 
55 counts in each of these phase bins. Using Poisson statistics, the 76 counts detected at phase 
0.35 deviate from the expected 55 counts with 99% confidence. Scattering and/or absorption of 
the x-ray emission could reduce the degree of modulation at either phase 0.15 or 0.35. Taking 
into account either possibility, the chance probability for the observed variation at phase 0.35 is 
2%. A decrease in the count rate is predicted at orbital phase 0.25 and the 44 counts we detect 
in this phase bin have a 4% chance probability of being drawn from a steady flux distribution. 
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